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A series of Yl^-^Sr^CruOs single crystals was measured to study the influence of site disorder 
on the transition line H SS (T) between quasi-ordered vortex lattice and highly disordered vortex glass, 
as well as on the maximum critical current density within the glass phase, j™ ax . When 32% of Ba 
is replaced by Sr, j" lax is an order of magnitude higher than in the unsubstituted compound. In 
contrast, the transition field H BB first drops by a factor of about five with a substitution of just 10% of 
Sr for Ba, and then remains approximately constant for higher Sr contents. Our results indicate that 
in very clean systems the order-disorder transition is affected very strongly by any crystallographic 
disorder, while above a certain threshold it is relatively robust with respect to additional disorder. 
In all substituted crystals H ss monotonically decreases with an increase of temperature. 

PACS numbers: 74.72.Bk, 74.25. Dw, 74.60.Ec, 74.60.Jg 



I. INTRODUCTION 

The vortex matter of cuprate superconductors has a 
rich phase diagram in the H-T plane (see e.g. Ref. [I]). 
Not long after the discovery of the cuprate supercon- 
ductors it has been established that there are at least 
two phases, a fluid (liquid or gas)a phase at high tem- 
peratures and fields and a solid (lattice or glass) phase 
at low temperatures and fields, ft was shown that the 
transition between these two phases is of Jjksf order, at 
least in low fieldsman fixL.|a_p£itical point.coQ Recently, 
it has been shownBLH3'QHEJIi3llj that there are two dis- 
tinct solid phases. In low fields the vortices, farm a quasi- 
ordered lattice in a so-called Bragg glass, liSEj which is 
stable against the formation of dislocations. In-higher 
fields they form a highly disordered, entangledtl solid. 
The nature of the high field phase is not well under- 
stood, ft could be anything from a highly viscous fluid 
to a vortex glass-with unbounded barriers against vor- 
tex movement .E3'E3 We will simply term it "glass" hence- 
forth. 

The high-field glass phase is caused by weak, random 
and uncorrelated disorder due to point-like defects, such 
as oxygen vacancies.EJ This type of disorder is present 
even in the most pure crystals of cuprate superconduc- 
tors. Strong, correlated disorder due to extended defects, 
such as twin boundaries or columnar irradiation defects, 
has got a different influence on the vortex matter and. 
can lead to additionaLohases, e.g. a "Bose glass" phaseta 
located at low nelds.EJ Correlated disorder may not be 
present in all crystals. In the following, when we sim- 
ply speak of "disorder" , we are referring to weak random 
point-like disorder. 

The pinning induced by the disorder has two effects on 
the vortex lattice.EHl Firstly, it localizes vortices by trap- 
ping them in low-lying metastable states, hardening the 



solid. Secondly, it promotes transverse wandering of vor- 
tices from their ideal lattice positions. The mechanism 
of the destruction of tJifi.JaJ±Lcfi_triggering the transi- 
tion to the glass phase, Ej'E2lEirE3c3 is very similar to the 
mechanism of melting at higher temperature. Similar to 
the way the lattice looses translational order and melts 
into a fluid, when thermal fluctuations become of the 
order of the lattice spacing, it looses translational or- 
der and "melts" into a glass, when the disorder-induced 
line wandering becomes of the order of the lattice spac- 
ing. In both cases the stability against the formation 
of topological defects, such as dislocations, is lost. The 
massive proliferation of dislocations accompanying the 
destruction of translational order leads to an entangled 
configuration of vortices, i.e. topological order is lost as 
well. In highly anisotropic superconductors, the lattice- 
glass transition rpaA£-|fall together with a decoupling be- 
tween the layers,B3'E3 of whicb_there is some evidence in 
Bi 2 Sr 2 CaCu 2 8+( 5 (BSCCO).Eil This is not expected. in. 
less anisotropic systems, like YBa2Cu30y_5 (Y123).e3'E3 
There is considerable evidence that the phase transition 
between lattice and glass, like the melting-lranailion be- 
tween lattice and fluid, is of first ordeaiSoo'O and is 
even a continuation of the melting lineE3 

Since the solid-solid transition is an order to disor- 
der transition with the disordered phase located at high 
fields, the transition line H SS (T) can be expected to shift 
to lower fields with the introduction of ad ditJo jj ia L-dis- 
order in the crystal structure. Calculationsc3EiloE3 of 
H ss , based on a Lindemann criterion, confirm this expec- 
tation, fndeed, there seems to be a correlation between 
the position of H ss and the purity of a sample, as can 
be seen by comparing different measurements on Y123 
crystals. Also, studies of the effect of electron irradiation 
on the transition between lattice and glass have found 
a systematic decrease of H ss with increasing irradiation 



doseJM 

An alternative way to introduce structural disorder 
into a system is to partially substitute an element of the 
compound. However, the generally varying level of oxy- 
gen within samples, which is difficult to measure and con- 
trol, constitutes a serious problem. YBa2Cu40s (Yf 24) 
is a compound very well suited for studying the depen- 
dence of H ss on chemically introduced disorder, because 
in this compound the oxygen content is fixed to 8 per 
unit cell. Another advantage of Y124 is that twinning, 
which leads to strong correlated disorder, does not exist 
in this compound.c3 Filially, Yf24 has got an interme- 
diate anisotropy 7 sw 12J23 which, as compared to Y123, 
shifts the glass phase down to fields more easily attain- 
able experimentally. To exclude or at least to minimize 
concurrent sources of disorder in the system it is advis- 
able to use single crystals for any measurements. 

If vortices are not ordered in a lattice, they can bet- 
ter adapt to the local pinning potential. Therefore, crit- 
ical current densities can be expected to be higher in 
the glass phase. High critical currents in high magnetic 
fields are desirable for many applications. Investigations 
of the disordered solid state and it's boundaries are thus 
very important in view of practical applications. Many 
irradiation studies found that critical current densities 
can increase by orders of magnitude upon irradiating a 
sample.E3 Studies on compounds with chemical disorder 
have also found increases of j c and even the development 
of a second peak.EJ 

In this work, we study the influence of a substitution 
of isovalent Sr for Ba in YBa2- a: Sr a ;Cu408 on the lat- 
tice to glass transition line H ss (T) , the irreversibility line 
Hi TT (T) and the critical current densities j c . In the pa- 
per, we are going to give evidence of an enhancement of 
the maximum critical current density in the glass phase 
j-max ky niore than an order of magnitude due to struc- 
tural disorder introduced by Sr substitution. The change 
of the vortex matter phase diagram due to the influence 
of relatively strong quenched disorder will be presented. 



II. EXPERIMENTAL 

Single crystals of Y124, with and without Sr substi- 
tution, were grown using a high-pressure technicuie, the 
details of crystal growth are reported elsewhereEJ Sin- 
gle crystals of YBa 2 - I Sr :E Cu408 with Sr content up to 
x = 0.64 or 32% were obtained. The crystals were 
checked with EDX and thek structure analyzed with Sin- 
gle Crystal x-ray analysis.EII The transition temperatures 
T c of the crystals were determined with a SQUID magne- 
tometer by temperature sweeps (both zero field and field 
cooled) in an applied field of 1 Oe. All samples chosen 
for further magnetic measurements had a well-resolved 
transition with a transition width (10% to 90%) smaller 
than 2K (see Fig. g). 

Measurements of the dc magnetization were performed 
with a Quantum Design MPMS5 SQUID magnetome- 
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FIG. 1: Normalized M(T) curves of the crystals used in this 
work showing sharp superconducting transitions. The curves 
shown were measured in a magnetic field of H = 1 Oe, under 
zero field cooled conditions. 



ter, additional measurements were performed with a non- 
commercial SQUID magnetometer with a sensitivity bet- 
ter than 10 -6 emu. All measurements were done with the 
magnetic field applied parallel to the c-axis of the crystal 
and generally the magnetization was measured at fixed 
temperatures as a function of the external field being 
swept up and down. After each change of the field, we 
waited for a time of 5 min before measuring the magnetic 
moment four times using a scan length of 2 cm. The 
relatively short scan length was chosen to ensure maxi- 
mum homogeneity of the applied field and was necessary 
to ensure a sufficient resolution of all features of inter- 
est discussed below. The wait time of 5 min was chosen 
to avoid effects of the initial relaxation of the field of 
the superconducting magnet without extending the time 
of measurements too much. Throughout the paper H 
denotes the applied field, i.e. no correction of the demag- 
netization was made. This is not a problem, however, 
since the magnetization is always much smaller than the 
external field except in the low field range, which is not 
the particular interest in this paper. 

The critical current density jjEi was calculated from 
the width of the hysteresis loop using an extended Bean 
modelH The irreversibility field H m was determined as 
the field where increasing and decreasing field branches 
of the M(H) hysteresis loop meet, or equivalently as the 
field where the critical current density vanishes. A uni- 
form criterion of j c = 50 Acm~ 2 was used. For the crys- 
tals studied it corresponds to a criterion of in between 
4.4 /iemu and 26 /xemu, which is not far from the accu- 
racy limit of our experimental setup. 



III. RESULTS AND DISCUSSION 

In the following we will discuss the influence of Sr sub- 
stitution on disorder and structure, critical current den- 
sities, the irreversibility line, and the transition from lat- 
tice to glass. 



A. Disorder and Structure 

Sr substitution increases the disorder in the structure, 
evidenced fo* example by a dramatically increased NQR 
line width. Ell No evidence for an introduction of corre- 
lated disorder was found and since twin boundaries are 
absent in YBa2- 2 ;Sr :r Cu408 the assumption of the dis- 
order in our crystals to be purely random and point-like 
seems justified. 

The Sr substitution also induces a variation of struc- 
tural parameters, causing a charge redistribution. Bond 
valence sum calculations and nuclear quadrupole reso- 
nance measurements indicate a transfer of holes from 
oxygen to copper atoms in the .G11O2 planes and from 
copper to oxygen in the chainsEil The charge redistri- 
bution can explain the increase of T c upon substituting 
more than 10% Sr for Ba. The initial-decrease of T c may 
be due to the increased site disorder EJ 

An important structural modification due to Sr substi- 
tution is a decrease of the thickness of the blocking layer 
db separating the superconducting CuC>2 planes, linear 
with Sr content. A substitution of 32% Sr_£pr Ba causes 
the blocking layer to shrink by about 1%,IEI which leads 
to an improved coupling between the superconducting 
planes. Torque measurements indeedpindicate a corre- 
sponding decrease of the anisotropy 7.EJ 



B. Critical Current Densities 

Figure |2p,) shows one of the half hysteresis loops mea- 
sured. The second peak is clearly discernible and dom- 
inates M(H). In all crystals, the second peak was ob- 
served over the whole temperature range measured, even 
just IK below T c . The form of the hysteresis loop and 
thus of j c {H) can be linked to the different phases of 
vortex matter. In low fields hysteresis is small and the 
hysteresis loop is rather anisotropic, which indicates that 
the major source of hyster esis in low fields is not of bulk 



origin (see also Sec. IIIC). The origin of the low bulk 



hysteresis is probably that the vortices are in the lattice 
phase in small H. Weak random point disorder is not 
effective in pinning an ordered dislocation-free lattice of 
vortices. The magnitude of the critical current density in 
low fields gives an indication of the degree of disorder in 
the system. For our samples, measurements of the critical 
current density at reduced temperatures of T/T c = 0.4 
and 0.7 in fields of 1 and 2kOe suggest that the disorder 
is increased monotonically up to the highest substitution 
level measured, with the biggest rate between x = 0.4 
and x — 0.64. 

When, upon rising the field, the lattice is destroyed in 
the glass state and the vortex system has an additional 
"dislocation degree of freedom" ,E2l individual vortices can 
much better adapt to local minima in the pinning po- 
tential. This results in a higher effective pinning force 
and thus a higher critical current density j c and relax- 
ation barrier .c3 Additionally, the entanglement itself was 
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FIG. 2: a) One of the hysteresis loops measured. The di- 
rection of the field change is shown by arrows. The onset of 
the second peak, the second peak maximum and the kink in 
M(H) between them are, denoted as H on , H2 P and H ss , re- 
spectively. Note that there are small differences between the 
respective fields between the two directions of field change 
(see text and Fig. hOJ) . The irreversibility field, where field in- 
creasing and decreasing branches meet, is denoted as Hi IT . b) 
Additional parts of hysteresis loops showing the kink region 
at three different temperatures. 



also suggested to increase critical current densities, due 
to an increased intervortex viscosityBj . prov ided the bar- 
riers to flux cutting are high enough.Eilca In our case, 
a steep rise of the magnetization occurs just below the 
kink field, denoted H ss in Fig. 0a) (see also Sec. IIID| ). 
Local magnetic measurements on a BSCCO crystal re- 
vealed a sharp increase inJA/(.ff)|, followed by a mono- 
tonic decrease afterwards!! The so-called "second peak" 
(or "fishtail") effect in BSCCO can thus be accounted 
for by the order-disorder transition. On untwinned Y123 
crystals, a similar sharp increase ipJAf (H)\ was observed 
by local magnetic measurementsjij but contrary to the 
situation for BSCCO, \M (H)\ continues to increase upon 
increasing field, although with a slower rate, which leads 
to a sharp kink in M(H). Only in even higher fields the 
magnetization starts to decrease with increasing field. 

Independently of vortex matter phase transitions, the 
second peak or fishtail effect in YI23 was attributed to 
a crossover from elastic to dislocation-mediated plastic 
creep. eIlH This mechanism is closely connected to the 
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FIG. 3: Critical current density at H^p vs. reduced tem- 
perature, for different Sr contents. It can be seen that un- 
substituted Y124 has much lower maximum critical current 
densities than substituted Y124, at all temperatures. Inset: 
jc{H 2p ) vs. Sr content x at T/T c = 0.7. 



lattice-glass transition though, since in the lattice phase 
dislocations are strongly suppressed while in the glass 
phase dislocations proliferate, as mentioned above.ES 
From collective creep theory a rise of the experimen- 
tally obsexjiable critical current density with the field is 
expectcdJilO as long as elastic creep can be assumed 
and the bulk pinning is relevant. From the shape of 
the hysteresis loops shown in Fig. [| follows that in 
YBa2_ a; Sr a ;Cu408 the situation is similar to the one in 
Y123. After rising the field even more, above H mi the 
vortex matter finally enters the fluid phase, where critical 
current densities are zero. 

The dependence of the hysteresis loop of an unsub- 
stituted Y124 single crystal on defects created by fast 
neutron irradiation was studied in Ref. [39| Fast neutron 
irradiation leads to improved pinning properties due to 
collision cascades (strongly pinning extended defects of 
spherical shape) and agglomerates of smaller defects. A 
clearly discernible second peak was observed in the un- 
treated crystal. In the weakly irradiated crystal the sec- 
ond peak was less pronounced and after irradiation to a 
fluence of 10 17 cm -2 it disappeared. At the same time, 
j c drastically increased in low to medium fields and the 
hysteresis loops became more symmetrical. 

The effect of fast neutron irradiation to remove the 
second peak may be understood from the different defect 
structure. The disorder resulting from the collision cas- 
cades is to some extent correlated and fast neutron irradi- 
ation can therefore create a Bose glass phase in low fields. 
Bulk pinning is also strong in the Bose glass phase, and 
if the irradiation induced correlated disorder is strong 
enough, the Bose glass phase may well destroy all of the 
Bragg glass phase. In high fields fast neutron irradiation 
is not as effective in enhancing j c and the irreversibility 
fields -ffi rr are not increased. 

In our case of Sr substitution the significant fields shift, 
but the qualitative form of the hysteresis loop remains the 
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FIG. 4: Dependence of the maximum pinning force density 
^p?n X = H2pjc(H2 P ) on Sr content, at a temperature of 65 K. 
The maximum pinning force density increases by more than 
an order of magnitude upon substituting 32% of Sr for Ba. 
Upper left inset: Second peak field vs. Sr content at T/T c = 
0.7. Lower right inset: Maximum pinning force density vs. Sr 
content, at T/T c = 0.7. 



same. As with fast neutron irradiation, Sr substitution 
increases j c , but most pronounced in medium fields, in 
the second peak region. The maximum critical current 
densities there, j™ ax (T) = j c {H = H 2pi T), increase dra- 
matically upon Sr substitution, as can be seen in Fig. 
g. At the temperature of 60 K, the critical current den- 
sity reaches almost the values of critical current densi- 
ties for Y123: For the crystal with 32% Sr substitution 
j c max (T = 60 K) = 6 x 10 4 AcuT 2 , while for a ^23 sin- 
gle crystal j c max (T = 60 K) = 6.3 x 10 4 Acm- 2 .El This is 
not the case, however, at the lower temperature of 50 K, 
where Y123 crystals have a maximum critical current 
density roughly three times higheiO than the one mea- 
sured on Sr substituted YBa2_ a; Sr 2 ;Cu40g. The main 
panel of Fig. || shows the dependence on Sr content of 
the maximum pinning force density F^ x — H2 P j c {H2p) 
corresponding to the critical current densities measured 
at T = 65 K. It can be seen that F^£* increases by 
more than an order of magnitude upon Sr substitution. 
The lower right inset of Fig. shows the Sr dependence 
of F^ x at a fixed reduced temperature. The initially 
approximately linear rise of the maximum pinning force 
density saturates for x rs 0.4 indicating that increasing 
the substitution level to x = 0.64 does not increase the 
pinning drastically any more. The critical current den- 
sity itself also shows saturation behaviour, but for lower 
Sr content (see inset of Fig. ||). The difference is due 
to the substitution dependence of the second peak field, 
shown in the upper left inset of Fig. 0. The dependence 
H2p(x) is similar to the substitution dependence of H ss , 
discussed in Sec. HIDj 



C. Irreversibility Fields 

The irreversibility fields H nl as a function of temper- 
ature are, for different Sr substitution levels, presented 
in Fig. |5j. The theoretical meaning of the irreversibility 
field is not very clear, often parts of it are controlled by 
surface or geometrical barriers (see below). In high fields, 
it was often found to be located near the transition from 
glass to fluid (see e.g. Refs. |l0|]l2] ). However, whether a 
true phase transition between a glass and a fluid phase 
even exists, is not a settled question. Not much is known 
about the properties of the highly disordered entangled 
solid located above -ff ss , apart from properties that dis- 
tinguish it from the lattice phase. The feature differenti- 
ating the phase from the liquid phase is a rapid freezing 
of the dynamics, but it is not sure whether this freezing 
of the dynamics is a proper phase transition or merely 
a crossover. The existence of a separate so-called vortex 
glass (VG) phase, distinguished by an unbounded dis- 
tribution of the heights of barriers between metastable 
states, was proposed early after the discovery of high T c 
superconductivity.tj Vortex glass scaling of current volt- 
age relations was used to determine the vortex glass tran- 
sition line H g experimentally, and also provided strong 
support for the VG theory (see e.g. Refs. [|[h]). How- 
ever, later experiments and numerical calculations found 
discrepancies (see Ref. O and references therein), and 
suggest a glass phase more akin to window glass. De- 
spite these uncertainties, H- m (T) is certainly important 
from a more practical point of view, as it limits the field 
range for applications. 

It is hard to see a clear trend in the Sr substitution 
dependence of the measured irreversibility fields. Differ- 
ences between H lrv of the different samples measured are 
relatively small, and the details depend on the temper- 
ature as well. Now it is important to remember that 
bulk pinning is not the only source of magnetic hys- 
teresis. There are additionaLeantributions due to Bean- 
LivingstpB surface barrierscilca and due to geometrical 
barriers O The significance of surface and geometrical 
barriers depends on several factors. One of them is the 
strength of the bulk pinning - the smaller the (bulk) crit- 
ical current density the bigger the relative importance of 
barriers. In YF^-xSr^C^Os the critical current densi- 
ties depend, as discussed above, strongly on x, and are 
very low for the unsubstituted crystal. Indeed, while 
the hysteresis loops of substituted crystals are generally 
rather symmetric, this is not the case for the unsubsti- 
tuted crystal, which is an indication that barriejr hystere- 
sis is more important than bulk hysteresis.c3 Another 
factor is the sample shape - bulk hysteresis depends lin- 
early on the sample diameter, while the barrier hysteresis 
scales with the ratio of thickness to diameter. Thus, bar- 
rier effects are more important for thicker and smaller 
(in diameter) crystals. The thickness to diameter ratio 
varies from 0.1 to 0.3 for the crystals used in this study. 
Finally, an important factor is the temperature, since 
the bulk hysteresis is depressed by thermal depinning, 




FIG. 5: The graph summarizes the dependence on the reduced 
temperature of the irreversibility fields (full symbols) and H ss 
(open symbols), for different Sr contents in Yf^-iSr^CXuOg: 
circles correspond to x — 0, diamonds to x = 0.2, down tri- 
angles to x — 0.35, up triangles to x — 0.4 and squares to 
x = 0.64. Error bars have been omitted for clarity. Full lines 
are fits of the H aa (T) data to Eq. (0). Except for x — 0, fit- 
ting was restricted to the range 0.65 < T/T c < 0.9. Dashed 
lines are guides for the eye indicating an exponential temper- 
ature dependence of H ss at lower temperatures. Variations of 
Hi IT (T) generally are relatively small and the details depend 
on the temperature. H sa (T) of the unsubstituted crystal is 
located in much higher fields than H SS (T) of all substituted 
crystals and for high temperatures is practically located in 
the region of the irreversibility fields. 



while the one due, to geometrical barriers is less affected 
by temperature.E-3 In the second peak region, the main 
part of the hysteresis is clearly of bulk origin, as fol- 
lows from the shape of the hysteresis loops (see Fig. 0) . 
However, since the total hysteresis observed is simply the 
sum of bulk and barrier hysteresis, the irreversibility field 
may still be determined by surface or geometrical barri- 
ers, since, the bulk hysteresis goes to faster at higher 
fields.EJO In high fields, bulk critical current densities 
were found to go to zero exponentially.!^ A linear rela- 
tionship between j c and ln(_ff ) is indeed followed in the 
crystals we measured. Figure k| shows this relationship 
for YBai.8Sro.2Cu408- The figure shows also, however, 
that there are systematic deviations below a certain tem- 
perature dependent threshold of critical current densities. 
For higher temperatures the deviation appears at larger 
critical current densities. This may indicate that the ir- 
reversibility lines are indeed influenced considerably by 
barrier hysteresis at higher temperatures. 

On the other hand, measurements of the influence of 
fast neutron irradiation on the high temperature part of 
the irreversibility field on a Y124 single crystal found a 
slight decrease of H m upon irradiating the crystal. EJ Ir- 
radiation may reduce Bean-Livingston type surface bar- 
riers, but those are, like bulk pinning, strongly temper- 
ature dependent and were generally not found to be im- 
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FIG. 6: Critical current density vs. logarithm of the magnetic 
field, for the sample with 10% Sr substitution. The relation- 
ship is linear for high critical currents. Below a temperature 
dependent threshold, deviations appear. Dashed lines extrap- 
olate the linear relationship to j c — 0. 



portant near T c . Geometrical barriers, on the other hand, 
should not change upon irradiation. As the irradiation 
did change the irreversibility field at high temperatures, 
it would indicate that -£/j rr at high temperatures can not 
be due to Bean-Livingston or geometrical barrier hystere- 
sis. 

However, as mentioned above, the importance of barri- 
ers strongly depends on the sample shape and we cannot 
be completely sure to what extent barriers influence iJj rr , 
especially in the case of the unsubstituted crystal, where 
bulk pinning is very weak. In order to explain the effect 
of the Sr substitution on the bulk properties, it is there- 
fore advisable to compare Hi„(T) at not too high reduced 
temperatures. A comparison of H irr (T) at T = 0.7 T c is 
shown in Fig. pi. The position of the irreversibility field is 
not changed much, but the general trend that i7i rr rises 
with x is consistent with the also rising maximum critical 
current density. The only irreversibility field that does 
not fit into this trend is the one of the unsubstituted 
crystal, which may be attributed to barrier effects that 
still are not negligible at this temperature for this weakly 
pinned compound. 

An increasing glass to liquid transition field can be 
justified theoretically in a similar way as the decrease 
of the lattice to glass transition field. Since the glass 
to liquid transition line depends on the relation between 
the pinning energy and the thermal energy, and because 
the pinning energy is increased upon introducing more 
disorder into the system, it would seem natural that the 
transition occurs at higher temperatures. 

The dependence of the glass to fluid transition on 
disorder was studied, for electron irradiated Y123, by 
Nishizaki et al.1^3 They measured both the glass tran- 
sition temperature T g resistively, and the irreversibility 
field -ffi rr , and found the two lines T g (H) and H irr (T) 
to coincide. Nishizaki et al. find the opposite disorder 
dependence of the glass to fluid transition, in their case 
H s and H[ Tr decrease with increasing fluence of electron 
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FIG. 7: Sr dependence of H ss at temperatures of 0.7 T c (open 
diamonds) and 0.4 T c (open circles), and of Hi IT (0.7T c ) (full 
diamonds). For both temperatures H aa drops significantly 
by a small substitution of Sr for Ba, but remains essentially 
constant upon further Sr substitution. 77i rr also drops by 
a small Sr substitution, but less significantly. It rises again 
upon further substitution and is higher for x — 0.64 than for 
x — 0. The x error bars have been omitted for clarity. 



irradiation. 

There are different possible explanations for the oppo- 
site behaviour of the glass to liquid transition line when 
the disorder is increased between the present measure- 
ments and those of Nishizaki et al.. 

Firstly, the disorder introduced by substituting 10% 
or more Sr for Ba in Y124 is much stronger than the 
disorder introduced by an irradiation of 2.5 MeV elec- 
trons with a fluence of up to 2 x 10 18 electrons/cm 2 in 
untwinned Y123. This can be seen by comparing the 
phase diagrams. The irradiated Y123 crystals show a 
first order melting transition in fields up to 50—100 kOe, 
while all of our Sr substituted Y124 crystals show no 
first order melting transition at all up to temperatures 
very near T c . Nishizaki et al. note that the decrease of 
H g upon increasing disorder they measured is consistent 
with a theoretical prediction made for systems with weak 
pinning.ca The disorder in our substituted crystals may 
be too big for this to be applicable. 

Secondly, the rise of i7i rr with Sr substitution could 
also be explained by a decreased anisotropy of the substi- 
tuted crystals, since H ll:r was found to scale with 7~ 2 s~ 1 , 
where s is the distance between two adjacent supercon- 
ducting planesB If the rise of H m in YBa2- a; Sr a; Cu408 
with x for x > 0.2 is due to a decreased anisotropy caused 
by the decreased blocking layer thickness, the initial drop 
for x < 0.2 may be due to the disorder, in which case the 
dependence of H m on disorder would be qualitatively the 
same as the one found by Nishizaki et al.. 



D. Transition from Lattice to Glass 



As noted in Sec. IIIB| , a kink in M{H) (denoted H ss ) is 
visible between onset and maximum of the second peak of 
the hysteresis loop shown in Fi g . M . A ccording to previ- 
ous studies performed on Y123,ll3ll3il3 it is likely that the 
kink corresponds to the transition line from the lattice to 
the glass state of vortices. In Ref. [13] it was argued that 
the equilibrium phase transition corresponds to the kink 
in M(H) in decreasing field, while the kink in increas- 
ing field corresponds to the upper limit of metastability 
of the lattice phase. However, since the pinning is much 
stronger in the glass phase, the equilibrium transition can 
be expected to be located quite close to the upper limit of 
metastability. This was indeed found from the measure- 
ments of Ref. [13| Also in our case, while we have more 
data on the field increasing branch, the kink in H de- 
creasing was checked at representative temperatures and 
the difference in the two kinks was found to be rather 
small indeed (see Fig. g). We should mention that the 
actual equilibrium transition is a little bit lower than the 
field of the kink in M(H^), which we call H ss (see Fig. 
no) . It is interesting to note that generally the kink posi- 
tions (for both directions of field changes) do not depend 
on the waiting time (we tested between 30 s and 10 rain) 
after changing the the field. Additionally the kink field 
location does not depend on the field change step size and 
not on the small inhomogeneities of the magnetic fields 
(~ 0.02%) tested by using different scan lengths, but the 
shape of the M(H) curve around the kink does depend 
on step size and scan length. On the other hand, the 
onset fields H on , defined as the fields where the absolute 
value of the magnetization has a minimum, have a slight 
tendency to shift to lower fields upon increasing waiting 
time, while the absolute value of the magnetization at 
the onset is decreased significantly upon increasing wait- 
ing time. Additionally the onset measured in decreasing 
fields can be much lower (up to 50% at low temperatures) 
than the onset measured in increasing field. Figure @ 
shows a selection of magnetization curves measured, in 
the region where the kink is located. The sharpness of 
the kink varies with temperature and depends on the 
sample as well. The inset of Fig. B shows the onset re- 
gion of one of the curves enlarged and the derivative of 
the magnetization (average of the four scans), which can 
help to locate H ss . Within a range of a few Kelvin be- 
low T c the unambiguous determination of H ss becomes 
difficult nonetheless. 

The increase of the absolute value of the magnetization 
around H ss observed in our measurements is hardly sharp 
enough to be called a "jump" . One reason for this is the 
spatial averaging of the inhomogeneous induction inside 
the sample. Additionally, the difference between the on- 
set of the second peak in increasing fields Hl n , defined 
as the minimum of the absolute value of the magnetiza- 
tion before the second peak, and the kink field H ss can 
be explained by the existencepof a region of metastabihty. 
around the phase transition.E3 Recent experimentstJ'Ej 
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FIG. 8: M(H) curves for increasing field near the onset of the 
second peak. Dashed lines are guides for the eye pinpointing 
the kink in the magnetization H SB . The inset shows a magnifi- 
cation of the onset region of a M(H) curve and it's derivative 
at T = 58 K. H ss is located at the position of the steepest 
increase of AM /AH . 



indicate that an abrupt change of the field injects a tran- 
sient disordered vortex phase at the sample edges. If 
the thermodynamically stable phase is the lattice phase, 
the transient disordered vortex phase then decays. This 
decay happens with a rate decreasing to as the field 
reaches H ss . In our measurements, the fields ilj n , H^ n 
and ifj-g was found to follow the same dependence on 
temperature and Sr content as H ss (HJ S ). Below the on- 
set, both the magnitude of the critical current density 
and the average value of the magnetization for the two 
branches of the hysteresis loop suggest that in this region 
surface and geometrical barriers are more important than 
bulk pinning in our samples, i.e. bulk pinning seems to 
be relevant only if the vortex matter is in the glass phase. 

The H SS (T) lines, as determined by the kink in the 
M(H') curves for the measured samples, are also shown 
in Fig. a. For all samples, H ss decreases monotonically 
with increasing temperature. A monotonically decreas- 
ing H ss was also found- .QiL-Mdi.ssCeo.isCuO^ap but 
measurements on Y123E£Ei}E3li3 showed an H ss increas- 
ing with temperature. The position of the phase transi- 
tion between quasi-ordered lattice and highly disordered 
glass is determined by the interplay between elastic and 
pinning energy. For the case of weak, random disorder 
due to point-like defects the thecpy of collective pinning 
was developed by Larkin et aLca The defects can in- 
teract with the vortices in two ways.EJ They can cause 
a spatial variation of the transition temperature (ST C - 
pinning), described by a modulation of the linear term 
of the Ginzburg-Landau free energy functional. Alterna- 
tively, they can cause a spatial modulation of the mean 
free path (<W-pinning) , described by a modulation of the 
gradient term of the free energy functional. In both cases, 
the influence of disorder is described by a disorder pa- 
rameter 7, proportional to the defect density. However, 
the temperature dependence of 7 is different for the two 
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cases. For <5T c -pinning, 7 ex 1/A 4 , while for ^-pinning, 
7 oc 1/(A£) 4 , where A and £ are the penetration depth 
and the coherence length£l The order-disorder transition 
position was c alcul ated analytically by using a Lindc- 
mann criterionB3lHi For the case of not too big anisotropy 
and disorder, following the calculation of Ref. 20, we get 



H SI 



H 



(1) 



with H Q = 2c 2 $ /£ 2 , U = $2 Cl £/(16V2tt 2 A 2 7) and the 
collective pinning energy U c = ((7$ 2 ^ 4 )/(167r 2 A 2 7 2 )) 1 / 3 . 
Here, $ is the flux quantum and cl ~ 0.1 — 0.2 the Lin- 
demann number. H ss is inversely proportional to both 
the anisotropy 7 and the disorder parameter 7. Con- 
cerning the temperature dependence, as long as we are 
below the depinning temperature Td p ~ C4/&B, we get 



H ss oc r 3 oc (1 - (T/T c ) 4 ) 3 / 2 



(2) 



in the case of <5T c -pinning and 



H ss oc £ ex (1 - (T/Te) 4 )- 1 / 2 (3) 

in the case of (W-pinning.E3 The calculation of Ref. |2l] 
leads to qualitatively the same results. Since H ss de- 
creases with temperature in our case, the pinning in 
YBa2_ E Sr a ,Cu408 cannot be S£ pinning. The temper- 
ature dependence of H ss of the unsubstituted crystal 
agrees satisfactorily with the formula proposed for 6T C 
pinning. In the case of the substituted crystals, how- 
ever, the agreement is limited to values of T/T c between 
about 0.65 and 0.9 (see Fig. ||, the full bold lines are fits 
to Eq. (g) ) . The lowered H ss in the vicinity of T c may be 
attributed to the finite transition width. 

However, at low temperatures, the temperature depen- 
dence is rather exponential in the case of the substituted 
crystals, as indicated by dashed lines in Fig. ||, with some 
indications of a flattening at the lowest temperatures 
measured. This behaviour is clearly at odds with Eq. 
(||) and to the best knowledge of the authors is not pre- 
dicted by any present theory of the order-disorder tran- 
sition. An exponential upturn in the onset field i?J n at 
low temperatures was also foundxm a Nd 2 _ x Ce2;Cu04_ x 
single crystal by Andrade et al.jt3 who attributed it to 
Bean-Livingston surface barriers. However, since the ob- 
served magnetic hysteresis is the sum of bulk and bar- 
rier hysteresis, it is difficult to imagine why the kink in 
M(H) at H ss should be influenced by surface or geomet- 
rical barriers, and in our case a clear kink is observable at 
low temperatures where H ss has already clearly departed 
from the (1 — (T/T c ) 4 ) 3 / 2 dependence observed at higher 
temperatures. 

A possible cause for a changed temperature depen- 
dence is a dimensional crossover. Due to the layered 
structure of the cuprate superconductors, vortex lines 
should be thought of being composed of stacks of "pan- 
cake" vortices.EJ Only below the 2D/3D crossover field 
H2D ~ &o/(j 2 s 2 ), where s is the distance between two 
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FIG. 9: Has and H\ TT in the vicinity of T c , for substitutions of 
10% and 32% Sr. 



adjacent layers, is the interlayer interaction between pan- 
cake vortices larger than their intralayer interaction and 
the pancakes form well defined vortex lines. For all our 
crystals, i?2D ^ 77kOe is estimated well above the up- 
per limit of the fields attainable in our magnetometers. 
Also, H ss cc £~5/ 2 (£5/ 2 ) for 2D and 5T C (5£) pinning,E3 
i.e. the temperature dependence should be even flatter 
in the 2D regime. A dimensional crossover can there- 
fore not be responsible for the observed upturn of H ss at 
low temperatures. In principle a second pinning mecha- 
nism, which is very effective at low temperatures, could 
be overlaid. However, point-like disorder should rather 
(additionally) suppress H ss , while the influence of corre- 
lated disorder should be more visible at high tempera- 
tures or low fields .Ej Figure g shows H ss and i?i rr of two 
of the substituted crystals in the vicinity of the transi- 
tion temperature. No sign of a tricritical point, where 
the two lines would meet, can be seen. A tricritical point 
would have to be located very near T c , where reliable 
measurements become increasingly difficult. 

It can be seen in Fig. || that any substitution of Sr low- 
ers H ss very significantly. For two values of T/T c , H ss vs. 
strontium content x is shown in Fig. |7J. After a substitu- 
tion of just 10% Sr for Ba, H ss drops roughly by a factor 
of 5, at both temperatures. A further increase of Sr sub- 
stitution, however, does not reduce the magnitude of H ss 
any more. Rather, H ss (x) remains essentially constant. 
We do not consider the peak at x w 0.4 as big enough to 
be significant. 

The initial decrease of H ss upon Sr substitution does 
not come unexpected. It follows from Eq. (|lj) that H ss oc 
(77) _1 . Sr substitution increases the d isorde r parameter 
and decreases the anisotropy (see Sec. |HI A| ). The large 
decrease of H ss upon substituting 10% of Sr indicates 
that the influence of the Sr substitution on disorder is 
much bigger than it's influence on the anisotropy. 

What is more surprising is the apparent saturation of 
the influence of the additional disorder for a substitution 
level x > 0.2, while the maximum pinning force is still 
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FIG. 10: Phase diagrams of the crystals with the lowest 
(x = 0, a)) and the highest (a; = 0.64, b)) Sr substitution. The 
quasi-ordered lattice phase at low temperatures and fields is 
separated from the highly disordered glass phase by the tran- 
sition line H aB (full down triangles, [thick line). Around the 
transition is a region of metastability ,E3 where both phases can 
coexist (shaded). The lower limit of metastability is marked 
by ffi„ (open down triangles) and the upper limit, very close 
to the equilibrium transition, by H BS (full up triangles). Also 
shown is Hj n (open up triangles). The irreversibility field 
Hi r r (full circles) is assumed to correspond roughly to the 
transition line between glass and fluid phases. Both H aa and 
Hi rl axe decreasing monotonically with increasing tempera- 
ture. For the crystal with x — the region of a possible 
tricritical point was not measured and is blurred in the fig- 
ure. For the crystal with x — 0.64 a tricritical point could not 
be detected and would have to be located very near T c . 



linearly increasing up to x ~ 0.4. In the case of weak 



point-like pinning H ss oc 7 



i dis , where 7 is the dis~ 



order parameter and ridis the defect density, is expected.El 
It would seem natural to assume the defect density to be 
proportional to the Sr substitution level x. The blocking 
layer thickness on the other hand, which is the structural 
factor with the biggest-influence on the anisotropy, de- 
creases linearly with x£2l A large initial decrease of H ss 
followed by an almost independence on x is therefore dif- 
ficult to explain. 

A downward shift of _ff ss upon introducing more dis- 
order was also found—in electron irradiated crystals of 
BSCCOQ and Y123.E3 The notable difference between 



both of these irradiation experiments and Sr substitu- 
tion in Y124 is that in both cases of electron irradiation 
there is no sign of any saturation of the disorder induced 
lowering of H ss . As discussed above, the disorder in- 
duced by Sr substitution is probably, even for the lowest 
substitution level of x ss 0.2, already much higher than 
the disorder induced by the electron irradiation with the 
largest fluence measured in Refs. and [L2| It may be 
that in our substituted crystals the disorder is already 
too big for the applicability of formulas derived for weak 
disorder. 

The resulting phase diagrams for YBa2Cu40s and 
YBai.3 6 Sr .64Cu4O8 are contrasted in Fig. 10. The 



glass phase of YBai.36Sr .64Cu4O8 has, as compared to 
YBa2Cu40g, expanded dramatically from a rather small 
part of the easily experimentally accessible phase dia- 
gram to the phase covering the biggest area in fields 
H < 55kOe. The most dramatic shift is the lowering 
of .Hss, which is to be expected as the (quasi-)ordered 
phase should become smaller upon introducing more dis- 
order into the system. The region of metastability is 
much smaller for the substituted crystal, indicating that 
strong disorder tends to reduce over-heating and espe- 
cially under-cooling effects. We speculate that this may 
be connected to our observation that for unsubstituted 
Y124 with very weak disorder H ss depends strongly on 
the exact amount of disorder, while for substituted Y124 
with rather strong disorder it almost does not depend on 
the exact level of disorder. In real crystals, the disorder 
density depends on the position on a mesoscopic scale. 
Soibel et alEn observed that for melting supercooling ex- 
ists only at local maxima of the transition field. In a first 
approximation the width of metastability can be linked 
to the difference between maximal and minimal transi- 
tion fields within the sample. In our case the same spatial 
variation of disorder should produce a much bigger vari- 
ation of the local transition field for the unsubstituted 
crystal than for the substituted ones. 



IV. CONCLUSIONS 

We investigated the influence of structural disorder, in- 
troduced by a substitution of Sr for Ba in YBa2Cu40s, 
on the critical current density at elevated fields and on 
the borders of the glass phase of vortices. We stress that 
the disorder introduced chemically by means of Sr sub- 
stitution is random and point-like, rather than disorder 
due to extended defects. This is shown by the hysteresis 
loops keeping their general form, similar to the case of 
electron irradiation and unlike the case of neutron irra- 
diation, where at least partially the effect is to introduce 
extended defects. 

Our main conclusion is that the introduction of even a 
small amount of disorder into very clean systems changes 
the phase diagram drastically by lowering the order- 
disorder transition line H SS (T) and also changing it's 
temperature dependence. When the disorder reaches a 
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certain threshold, however, introducing additional disor- 
der does not continue this tendency. With other words, 
in highly disordered systems (including probably most 
systems with partial chemical substitutions), the vortex 
matter phase diagram is relatively robust with respect to 
variations in the exact degree of disorder. 

This conclusion is supported by the rather small vari- 
ation of the measured kink location as an indication of 
the order-disorder transition at fixed reduced tempera- 
ture for any Sr substitution level x > 0.2, especially as 
compared to the corresponding transition in a clean, un- 
substituted crystal. That the disorder is still increasing 
with increasing substitution level for x phO.2 is indicated 
both by an increasing NQR line widthEil and increasing 
critical current densities at low fields and the maximum 
pinning force density. 

In strongly disordered systems, metastability is less im- 
portant than in very clean systems and especially the re- 
gion of undercooling of the glass phase is much smaller. 
This might be linked to the observation that the position 
of the transition line does not depend much on the exact 
amount of disorder in strongly disordered systems, since 
variations of the disorder density on a mesoscopic scale 
do not produce large positional differences of the local 
transition field in this case. 

Similar to the robustness of the phase diagram above a 
disorder threshold, both the maximum pinning force den- 
sity and the maximum critical current density in the glass 
phase increase initially with increasing disorder, the crit- 
ical current density increases by an order of magnitude 
upon substituting 10% Sr. But they then show satura- 
tion behaviour, i.e. the pinning force density and critical 
current densities can be raised only up to a certain point 
by introducing random point-like disorder. This may be 
relevant for possible applications where high critical cur- 



rent densities are required. 

It should be noted that current detailed theories of the 
order-disorder transition generally assume the disorder 
to be very weak. Our measurements indicate that in sys- 
tems with a high amount of disorder, the phase diagram 
and in particular the order-disorder transition differ qual- 
itatively from the weak disorder case. Theoretical inves- 
tigations beyond the weak disorder limit would be very 
helpful. 

Both the pinning of unsubstituted Y124 and Sr sub- 
stituted Y124 is likely to be of 6T C , rather than St type, 
indicated by the temperature dependence of the order- 
disorder transition at intermediate temperatures. For 
the substituted crystals this can be expected, since there 
are variations of T c upon Sr substitution. However, the 
exponential-like temperature dependence of the order- 
disorder transition of substituted crystals at low tem- 
peratures is also emphasized. The dependence cannot be 
linked to the influence of surface barriers, at least not 
in any straightforward way, the same dependence of all 
significant fields (onsets and clearly discernible kinks for 
both field directions) suggests it is a true bulk property of 
the transition. We stress that the observed temperature 
dependence is at odds with currently published theoret- 
ical formulas. 
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